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Abstract 

Wearless friction is investigated in a simple mechanical model called Frenkel- 
Kontorova-Tomlinson model. We have introduced this model in [Phys. Rev. 
B, 53, 7539 (1996)] where the static friction has already been considered. 
Here the model is treated for constant sliding speed. The motion of the 
internal degrees of freedom is regular for small sliding velocities or weak in- 
teraction between the sliding surfaces. The regular motion for large velocities 
is strongly determined by normal and super harmonic resonance of phonons 
excited by the so-called "washboard wave" . The kinetic friction has maxima 
near these resonances. For increasing interaction strength the regular mo- 
tion becomes unstable due to parametric resonance leading to quasistatic and 
chaotic motion. For sliding velocities beyond first-order parametric resonance 
bistability occurs between the strongly chaotic motion (fluid sliding state), 
where friction is large and a regular motion (solid sliding state), where fric- 
tion is weak. The fluid sliding state is mainly determined by the density of 
decay channels of m washboard waves into n phonons. This density describes 
qualitatively the effectiveness of the energy transfer from the uniform sliding 
motion into the microscopic, irregular motion of the degrees of freedom at the 
sliding interface. For a narrow interval of the sliding velocities we also found 
enhanced friction due to coherent motion. In the regime of coherent motion 



non-destructive interactions of dark envelope solitons occur. 
PACS numbers: 46.30.Pa, 05.70.Ln, 81.40.Pq 
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I. INTRODUCTION 



The sliding of two solid bodies relative to one another is a non-equilibrium process where 
the kinetic energy of the uniform motion is transferred into energy of irregular microscopic 
motion, i.e. into heat. This dissipative process leads to friction, called dry friction, which 
differs from viscous friction. Since more than 200 years the phenomenological laws of dry 
friction (Coulomb's laws) are well-known and are well established in applied physics [|TJ. 
They state that the friction force is given by a material parameter (friction coefficient) 
times the normal force. The coefficient of static friction (i.e., the force necessary to start 
sliding) is larger or equal to the coefficient of the kinetic friction (i.e., the force necessary 
for sliding at a constant velocity). 

One cannot expect that Coulomb's law can be derived as rigorously as the laws of viscous 
friction. The reason is that the latter one can be deduced from the well-established theory 
of near-equilibrium thermodynamics whereas dry friction works mostly far from thermal 
equilibrium. In fact Coulomb's laws oversimplify a very complex behavior which involves 
elastomechanical, plastic, and chemical processes operating on different length and time 
scales (for a general overview of the interdisciplinary aspects of dry friction see Ref. |2j). 
Experimentally deviations of Coulomb's law are often found 0. For example, the kinetic 
friction depends on the sliding velocity. A common pattern is that it first decreases, then 
goes through a minimum and finally increases |§. In the case of a thin lubrication layer 
(thickness: a few monolayers) one often finds the contrary: The kinetic friction first increases, 
then goes through a maximum and finally decreases @]. 

In the last ten years the interest in dry friction (i.e., friction between solids) has been 
aroused again due to the possibility of reproducible friction experiments on the mesoscopic 
and nanoscopic scale (nanotribology) 0. There is some hope that for such less complex 
systems an understanding is possible by theoretical modeling ||]. From the general point of 
view of nonequilibrium thermodynamics there is also some hope that it should be possible to 
develop a simple model which shows all the basic phenomenological features of dry friction. 
Like the Ising model in the theory of ferromagnetism it is acceptable that such a model is not 
realistic in every details. What is more important is that it is on the one hand sufficiently 



complex to be a testing ground of basic concepts of the statistical mechanics of dry friction. 
On the other hand it should be simple enough to get analytical results. 

The Frenkel-Kontorova-Tomlinson (FKT) model goes back to an old idea of Tomlinson 
which explains Coulomb's laws on the atomic level |||7|. His attempt concerns the most 
puzzling feature of dry friction, namely, that the friction force remains finite in the limit of 
quasistatic sliding (i.e., sliding velocity going to zero). It is clear that this behavior cannot 
be treated, like viscous friction, as the linear response of a system near thermal equilibrium. 
Tomlinson's idea was that friction is caused by mutual "plucking" of surface atoms. Each 
plucked surface atom vibrates. Plucking is possible only if the sliding brings surface atoms 
into metastable states from which they suddenly jump into more stable ones. The energy 
difference is turned into kinetic energy which is totally dissipated due to the emission of 
some type of bulk or surface waves (e.g., sound waves) by the plucked atom. This is the 
basic dissipation channel. Tomlinson's mechanism works even in the limit of quasistatic 
sliding provided that the barrier between the metastable and the stable states is sufficiently 
high and the temperature sufficiently low. 

This simple picture assumes (i) that the relaxation time of a plucked atom is much smaller 
than the averaged time between two jumps of the same atom and (ii) that a vibrating atom 
does not excite vibrations of other surface atoms. These assumptions lead to a kinetic 
friction that does not depend on the sliding velocity ||||. 

The first assumption does not hold if either the vibration after the jump is not totally 
dissipated until the next plucking event of the same atom (weak dissipation) or the jump itself 
occurs on a time scale larger than the averaged time between two plucking events of a single 
atom (strong dissipation). One can imagine that in the case of atomically flat surfaces such 
plucking events occur in a regular fashion with a frequency (so-called washboard frequency |5|) 
which is determined by the sliding velocity multiplied by the principal wave number of the 
surface structure. In the case of weak damping we can therefore expect resonances which 
should lead to a velocity-dependent kinetic friction with resonance peaks. The situation 
can be modeled by the independent oscillator or Tomlinson model |]7|JlO|. It is a damped 
harmonic oscillator (eigenfrequency uq) interacting with a sliding surface (velocity v) which 
is described by a spatially periodic potential with wave number k. Two types of resonances 



are possible in the Tomlinson model JTTJ: Normal resonance where the oscillator is plucked 
after n oscillations (resonance condition: uj = vk = u /n), and parametric resonance which 
is an instability where oscillations are excited due to the modulation of the effective frequency 
of the oscillator during sliding (resonance condition: uj = 2uo/n, where n is an integer) ||12|| . 
In the case of weak interaction these resonances can be calculated analytically, whereas in 
the case of strong interaction the motion of the oscillator becomes chaotic and pronounced 
hysteresis occurs [10,|lT|. The kinetic friction becomes a wild, multi- valued function of the 
sliding velocity. This is even the case after averaging over an ensemble of identical oscillators 
0- 

The easiest way of taking into account the coupling between surface atoms is to introduce 
a linear nearest-neighbor interaction of the oscillators. This leads to the Frenkel-Kontorova- 
Tomlinson model which we have already introduced in Ref. ||. The FKT model assumes two 
rigid sliding bodies where the layer of surface atoms of one body is modeled by a mechanical 
particle-spring model (see Fig. |1]). There are two kind of springs, a nearest-neighbor spring 
and a spring which couples each particle to the body. The interaction between the two 
bodies is described by a spatially periodic potential. The dissipative coupling of the layer 
with the environment is modeled by damping terms which are proportional to the relative 
velocities. The static properties (ground state, meta-stable states, static friction) of the one- 
dimensional FKT model have already been discussed in Ref. ||. In this paper we investigate 
dynamical properties such as the kinetic friction of the one-dimensional FKT model. A brief 



summary of Ref. M and of this paper is given in Ref. [14 



The FKT model is far from being realistic in certain aspects. Its main disadvantages are 
the following. The assumption of harmonic particle interaction does not hold in the case of 
strong distortions. The elasticity of the bulk is missing which leads to an artificial gap in 
the phonon dispersion relation (see Sec. |1V A]) . 



The FKT model shares many similarities with the Burridge-Knopoff (BK) model IJT5 



which was proposed as a model of an earthquake fault. There also a linear, harmonic chain 
of blocks attached via springs at a rigid plate is assumed. But the interaction with the other 
surface is modeled by phenomenological force describing a dry friction law. Mostly velocity- 
weakening laws are used. That is, the friction force decreases with the sliding velocity. 



The dynamical properties of the BK model have been studied extensively in the literature 

EMI- 

The paper is organized as follows. In Sec. the FKT model and the kinetic friction 
Fx are defined. The behavior for small sliding velocities is treated in Sec. |m[ Here mostly 
stick-slip motion of individual particles occurs. At higher velocities excitations of phonons 



become important. The consequence of several kinds of resonances is discussed in Sec. |TV 



In the last section we compare our results with other models, especially with the BK model. 



II. DEFINITIONS 

The Frenkel-Kontorova-Tomlinson (FKT) model has already be introduced in Ref. ||. 
Thus the key features are described only briefly. Here we will concentrate more on the 
dynamical properties of the model which have not been investigated in Ref. 0. Also a 
definition of kinetic friction Fx is given. 



A. The FKT model 

A linear chain of N particles coupled harmonically forms the backbone of the FKT model 
(see Fig. [I]). The particles can move only parallel to the sliding direction. Thus the model 
is completely one-dimensional. All particles have the same mass which defines the mass 
unit. Also the nearest-neighbor interaction is assumed to be uniform which defines the force 
unit. The chain is the interface between two rigid sliding bodies which are modeled in two 
very different ways. The interaction with the lower body is described by a spatially periodic 
potential. Its periodicity defines the unit of the length scale. The interaction with the upper 
body is assumed to be harmonic. That is, each particle is coupled via a spring to the upper 
body. All springs are identical and they are attached with equal spacing at the upper body. 
The boundary conditions are assumed to be periodic. 

In this paper the motion of the chain is investigated for the case of sliding the upper 
body against the lower one at constant velocity v. This externally constrained sliding leads 
to vibrations within the chain. These vibrations excite waves into the bulk of the sliding 
bodies which carries energy away from the sliding interface. In the FKT model the degrees 



of freedom which are responsible for these waves are not explicitly present. We include 
their effect by introducing two phenomenological damping forces which are proportional to 
the particle velocity measured relative to the lower and the upper body, respectively. In 
this paper we investigate the FKT model in the limit of zero temperature where fluctuation 
forces disappear. Below we will see that the kinetic friction is directly given by the rate of 
energy which is transferred into these basic dissipation channels. 
The equation of motion of the FKT model reads 



where £j is the position of particle j relative to its equilibrium value given by the upper 
body, c is ratio of the lattice constant of the upper body to that of the lower body, b and 
k are the strengths of interaction with the lower and upper body, respectively, and ji and 
7(7 are the damping constants of the particle motion relative to the lower and upper body, 
respectively. The boundary condition is 
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which requires 



c = 



— , with M integer. 
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The transformation 
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turns (HD into 
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with 



7 = 1l + 7a- 



(6) 



In addition to the symmetries of the FKT model already mentioned in Ref. 0, eq. (|]) 
is invariant under translation in time by an integer multiple of the inverse sliding velocity. 
Below we will see that all solutions which do not break this symmetry are defined by a 
single one-parameter function, the dynamical hull function, which is similar to the static 
hull function 0. 

For the results presented in the following sections, the equation of motion fl5|) is often 
integrated numerically We have done this by using a second-order Runge-Kutta scheme 
which in the non-damping case is identical with the Verlet algorithm p5 |. For the time 
step we have chosen 0.01. This is sufficient because all simulations are done for k — 1 and 
v < 1.5. The corresponding fastest time scale is 2^/^/5. 



B. The kinetic friction 

The kinetic friction Fk is the temporal average of the force necessary to keep the body 
sliding at a constant velocity v. It is given by the sum of the forces between the particles 
and the upper body f5|, 

rp N 

F K = -\im^ jf £&(f)dt. (7) 

There is a second equivalent definition which clearly shows the dissipative nature of the 
kinetic friction. The idea is that after the decay of transients the energy flows through the 
system at a constant rate. That is, the energy which is put into the system per unit time 
by sliding is on average totally dissipated. Thus the total energy (i.e., the kinetic and the 
potential energy) of the chain is on average constant. With the further assumption that the 
motion of the £'s is bounded we get 

F K = lL vN + hm -L jT £g (t) dt (8) 

From the last definition we clearly see that the kinetic friction is proportional to the rate 
of energy transfer from the motion of the interface particles into some types of bulk waves. 
This is a general principle: The kinetic friction is determined by the various rates of energy 
transfer starting from the uniform macroscopic motion of the sliding bodies and ending in 
some dissipation channels. For the FKT model we are mainly interested in the efficiency of 



the transfer from the macroscopic sliding into microscopic motion of the interface particles. 
In section |IV] we will see that this efficiency is rather high if phonons are excited resonantly. 
That is, the kinetic friction as a function of the sliding velocity shows resonance peaks. 

In the following sections we always discuss the kinetic friction without the first term of 
definition (g). That is, in any plot of F K versus v one should add ^l^N . We will see that 
for very large velocities this term dominates the second one. 

III. THE LOW- VELOCITY REGIME 

In this section we investigate the FKT model for very small sliding velocities. In fact 
we are treating the quasistatic limit v — > 0. In Ref. f| we have already discussed the 
disappearance of the kinetic friction in this limit if b is less than some critical value . Here 
we investigate F K as a function of v for b below and above this threshold. 

A. b < b? 

In Ref. || we have shown that for b < b^ the stationary values of the £'s follow the 
quasistatic sliding adiabatically. Furthermore, there is only one stationary state which can 
be described by the static hull function g$ [defined by Eq. (10) in Ref. || which is identical 
with (|37D for v = 0}. Therefore, the solution of the equation of motion (|j) in the quasistatic 
limit reads 

= 9s(vt + cj) with g s (x + 1) = g 8 . (9) 
Using (^) we immediately see that the kinetic friction behaves viciousely, i.e., 

F K = JeftVN, (10) 

with an effective damping constant 

7« ff = 7x + 7 / (^) dx. (11) 

In the Tomlinson limit (i.e., k, b — > oo) and for c = 1 an analytic expression can be obtained 
(see appendix |A|). In the limit b — » b^ the effective damping constant diverges like 

7efr ~ (6f - b)- 1 ' 2 . (12) 



B. b > b* 



For b > additional metastable states appear P5||. Due to quasistatic sliding these 



states eventually annihilate with an unstable counterpart in a saddle-node bifurcation ||. 
At those points the chain has to find a new metastable configuration. That is, the particles 
of the chain no longer move adiabatically. They start to move with finite velocities, and 
after a transient they reach a new metastable state. The dissipated energy is the difference 
in potential energy before and after the reconfiguration event. This is the generalization of 
Tomlinson's friction mechanism. The kinetic friction in the quasistatic limit is therefore the 
averaged energy drop per jump (AV) divided by the averaged distance (Ax B ) between the 
relative positions of the sliding bodies between two subsequent reconfiguration events: 

F *<"-°> = ©- (13) 

In the overdamped limit, i.e., 7 — * 00, the system selects via steepest descent the nearest 
possible metastable state on the energy surface. Starting with an arbitrary metastable 
state the system will reach after several jumps a state which can be obtained by adiabatic 
deformation from the ground state of the undriven FKT model [J27J. Then this state always 
reappears after each jump. The reason for this selection rule is closely related to the fact that 
the ground state has the highest depinning force of all metastable states ||. The sequence 
of reconfiguration events becomes therefore regular. We assume that c is the ratio of two 
coprime integers P and Q (i.e., c = P/Q) where Q is a divisor of N. In each reconfiguration 
event N/Q equally spaced particles jump into the neighboring potential wells. Each jumping 
particle together with the neighboring particles dissipates AVo. Thus the totally dissipated 
energy is AVqN/Q. The distance between two events is l/Q. Thus 

F K (v 0) = AV N. (14) 

Near the onset of friction F% scales like 

(ir,) 



N 

The leading orders of this expansion are calculated in Appendix [A] for c = 1 and for the 
Tomlinson limit. For & > k > 1 we get (see appendix |A|) 



^-M^(?) § +-)- 

In the limit n/b — > the kinetic friction is approaching the static friction F$ = F™ 11 * = bN. 

Figure |2| shows for several values of c the static friction F$ and the kinetic friction 
Fk{v — > 0) which are obtained numerically by solving the delay equation of the hull function 
gs- First of all, we notice the excellent agreement between analytical and numerical results 
for c = 1. Furthermore, we observe that the kinetic friction is approaching the static friction 
from below not only for b — > oo but also for Q — > oo (i.e., for decreasing commensurability of 
c). In fact the difference scales like 1/Q for Q ^> 1. This can be understood by the following 
argument. The jumps of the N/Q particles represent a kind of microslip leading to typical 
stick-slip motion with a sawtooth curve of the lateral force F = —k J2j £j as a function of the 
relative position between the sliding bodies. The maximum of this curve defines F$ whereas 
the kinetic friction is the average in accordance with definition (0). Thus F$ — Fk(v — ► 0) 
is given by half the height of the jumps which is proportional to 1/Q because in each jump 
N/Q particles are jumping over a distance which is independent of Q for Q ^> 1. The jump 
is roughly proportional to the largest gap in the hull function gs (see Figs. 2 in Ref. 0). 

If we take inertia into account, i.e., 7 finite, the system does not necessarily select the 
nearest possible minima of the energy landscape. That is, the particle may jump over several 
potential wells before being captured. This has two consequences: (i) It lowers the kinetic 
friction which can become zero in the extreme case of 7 — > 0. (ii) It makes the motion 
irregular. That is, the states between the reconfiguration events are no longer equivalent to 
the ground state, and AV and Axb are fluctuating properties. 



This irregular behavior is similar to the irregular motion in the BK model |L7[. For 
example, as in the BK model the FKT model shows a power law distribution of AV up 
to a certain cutoff value ||28|| . This cutoff monotonically decreases with increasing k. For 



k ^ 1 the power law disappears and the distribution is exponential p8| , |24|| . The power 
law behavior has also fed the speculation that dry friction may be a self-organized critical 



process. This subject is not addressed here. The interested reader is referred to Ref. f28 
where the FKT model is discussed in this context. 

The dynamical behavior of the chain does not change qualitatively if the sliding velocity is 



finite but still small. There is again a regular sequence of reconfiguration events where N/ Q 
equally spaced particles jump. It is expected that F^(v) increases with v because the critical 
slowing down of the system near a saddle- node bifurcation delays the reconfiguration event. 
Thus the energy drop AV is larger than in the quasistatic limit because at the beginning 
of the reconfiguration event the strings connecting the particles with the upper body are 
stretched more. From simulations we find that dFx/dv |„^o is indeed positive (see Figs. ||| 
and ||). We also see that after reaching a maximum Fk decreases. This can be understood in 
the following way. In the simulations shown in Figs. § and |^ the motion is not overdamped. 
Therefore the particles oscillate after jumping. If the sliding is sufficiently fast they are still 
oscillating at the next jumping event. Thus they explore the phase space and may therefore 
overcome the barrier before it disappears in the saddle-node bifurcation mentioned above. 
For increasing sliding velocity the delay caused by critical slowing down is more and more 
compensated by this mechanism. It leads to a local maximum of Fk{v) at a velocity which 
scales like 7. This is roughly confirmed by the friction curves in Fig. [| 

IV. THE HIGH VELOCITY REGIME 

In the previous section we have investigated the FKT model in a regime where Tom- 
linson's basic dissipation mechanism works, namely, that the energy put into the system is 
released and subsequently dissipated in recognizable reconfiguration events. For increasing 
sliding velocity the average time between successive reconfiguration events will eventually 
become smaller than the average time scale on which such events take place. Reconfiguration 
events overlap and are no longer recognizable. Thus we have entered a new regime where 
the particles are always in motion. In this regime the friction may be strongly enhanced due 
to resonances of waves, i.e., phonons. 

In this section the consequences of these resonances are investigated on various levels. 
First we treat the linearized equation of motion. In a second step we calculate numerically 
all regular solutions which do not break the discrete temporal translation symmetry of the 
equation of motion. Then irregular solutions are investigated. At the end of this section we 
report on solitary behavior. 



A. The basic resonances 



The external potential describing the hard surface of the lower body makes the equation 
of motion (^) nonlinear. Therefore it is very hard to find analytical solutions. For a weak 
potential, i.e., feci, and sufficiently strong dissipation one can expect that the oscillations 
of the particles are small. Thus we can linearize the equation of motion, 

ij+iij = 0+1 -(2+«)& + 

+6sin27r(ft + cj) 

+27r6cos27r(wt + cj) ■ (17) 

Both terms on the right-hand side which are proportional to b can be considered as driving 
terms. These driving terms are periodic in t and j with frequency 27rv and wave number 
2nc. We call this driving wave the "washboard wave". It leads to different types of reso- 
nances. The first driving term is responsible for the main resonance. The last term leads 
to parametric resonance. Both terms are responsible for superharmonic resonance. 

The washboard wave can resonantely excite phonons (i.e., sound waves) which are solu- 
tions of ([D]) for 7 = 6 = 0. They have the form 

ilkmj±u(km)t] 
e J 

with the dispersion relation 



u(k) = Jn + Asm 2 ^. (18) 



Because of periodic boundaries, k has to be discrete 



2-7T 

k = —771, 771 = 0, . . . , JV - 1. (19) 



1. The main resonance 



We drop the last term of (|T7|). Such linearized equations have already been studied in 
many other models ]2"9"|-pT|] . They can be easily solved by the Ansatz 

= A cos[2n(vt + cj) + 0] (20) 



which leads to a typical resonance line 



A= . (21) 

y/[u) 2 (27Tc) - (2tto) 2 ] 2 + (27TV 7 ) 2 



Note that (pO) is the solution in the limit t — > oo where all other solutions have died out 
because of 7 > 0. 

In order to calculate from this solution the kinetic friction one has to be aware that for the 
linearized equation of motion ( fTTD the two definitions (|7|) and (|8|) are no longer equivalent. 
The first definition ( PUD yields zero friction whereas the second definition leads to 

Fk{v) = 27TV& 2 

N [uj 2 (2nc)-(2TTv) 2 ] 2 + (27rv 1 ) 2 ' { ' 

This result is quadratic in b because (§) is quadratic in £. To get the same result with (|7|) 
one has to calculate £ 3 - up to second order in 6 (see appendix ||). 



For u — > the kinetic friction ( 22"D goes to zero. Thus we have the frictionless case 
b < bf from Sec. A| . For c = 1 we recover in leading order of b the exact result (|X3|). The 
solution ( p2[ ) increases with v, reaching the maximum 

f = 4 + < 23 > 

at 

u = v? = + 0( 7 2 ) (24) 

and decreases to zero like 1/v 3 for v — ► 00. The resonance peak occurs at that washboard 
frequency 27rt; which corresponds to the frequency of the phonon having the same wave 
vector as the washboard wave. 

Qualitatively the same behavior of the kinetic friction as a function of the sliding velocity 
is often found in dry friction experiments between atomically flat surfaces with a few mono- 
layers of lubricating molecules in between ||. Usually such friction curves are interpreted as 
a dynamic phase transition from liquid-like behavior at small velocities to a solid-like state 
at high velocities. At the maximum of the kinetic friction the lubrication layer is assumed 
to be in an amorphous state P,|33[ . 

We can give an estimate of the validity of the kinetic friction (p2"D. The linearized equation 
of motion works fine as long as 1. Using ( pl| ) we find 



b 2 < [u 2 (2nc) - (2wvyY + (2vrw7) 2 . (25) 



The approximation (|22j) first fails at the resonance peak. Thus (|23|) is a good approximation 
as long as b -C u;(27rc)7. For sliding velocities much larger than y/b/fiir) the approximation 
is always good because the particles make only tiny oscillations around their equilibrium 
positions. We call this state the solid-sliding state. 

2. Superharmonic resonances 



The last term in (|T7|) is responsible for resonances of higher harmonics. For example, 



the solution (p0|) would lead to a term of the form 

^4ni(vt+cj) 

It is therefore natural to expand £j into a Fourier series: 

oo 

0= £ A m e 2 ™ m ^\ A_ m = A* m . (26) 

m=— oo 

The linearized equation of motion turns into a set of linear algebraic equations for the Fourier 
coefficients A m : 

ib 

LmA m — 7ib(A m _i + A m+ i) = — (5 m -i — 5 mj i) , (27) 
where 5i j is the Kronecker symbol and 

L m = uj 2 (2Trcm) — (2TTvm) 2 + 2irvm'yi. (28) 



For small values of b, ( P7|) can be solved perturbatively (see appendix |B|). In order to express 
Fk in terms of the Fourier components A m , we use again definition (||]) of the kinetic friction 



= 8vrV \^A m \ 2 . (29) 

^ V m=l 

Because A m is proportional to 1/L m we expect a peak in Fk(v) near the superharmonic 
resonance velocity 



This justifies an approximation of Fk which takes only the leading term in the power series 
of A m into account: 

2 



« 2 V1 Yi m2 n 

JV m=l j=l 



7r6 



(31) 



This is not a correct expansion of Fk in powers of b but it reproduces the superharmonic 
resonance peaks qualitatively very well as is shown in figure |](a). It also gives the right 
scaling of the height of the resonance peaks, namely, b 2m for the superharmonic resonance 
of order m. 

A special case of superharmonic resonance is m = 1 which is the main resonances. It 
can be shown that 

vf > v^, for m > 1. (32) 

Thus the superharmonic resonance peaks occur between zero and vf. But they do not order 
sequentially (i.e., ff < for some n > m) even though — > for m — > oo. 



3. Parametric resonances 

The last term in the linearized equation of motion (O) is also responsible for parametric 
resonance of phonons. Parametric resonance is an instability phenomenon. The prototype 
is a pendulum with a vertically oscillating support The linearized equation of motion 
is the well-known Mathieu equation [55| which is a (damped) oscillator with eigenfrequency 



u . Energy can be put into the oscillator due to periodic modulation (with frequency uS) of 
a parameter which determines u Q . If the oscillator is damped this energy will dissipate. 

The behavior of such a parametrically driven oscillator is very different from an oscillator 
which is driven by an additive force. In the latter case the response is linear whereas in 
the first case it is strongly nonlinear. The linear relationship between the amplitude of the 
oscillation and the strength of an additive force is caused by the fact that the driving force is 
independent of the amplitude of the oscillation whereas in a parametrically driven oscillator 
the force is proportional to the amplitude. This leads to a positive feedback mechanism. For 
a small proportionality constant (i.e., strength of driving) more energy is dissipated than is 
gained from driving and the oscillation dies out. But if the strength of driving exceeds a 



certain threshold, more energy is gained than is dissipated and the amplitude increases. In 
a linear equation like the Mathieu equation such solutions go to infinity. But in real systems 
like the pendulum, nonlinearities will eventually saturate this growth. 

The threshold for parametric resonance is a function on the frequency ratio ujq/uj. It has 
relative minima near the parametric resonance condition 

uj = — -, m=l,2, ... (33) 
m 

where m denotes the order. The main difference between a parametrically driven oscillator 
and the FKT model is that in the latter case there are several oscillators (i.e. phonons), 
and the driving is spatially non-uniform because it is a wave. In this case the parametric 
resonance condition reads (see appendix ||) 

v P {q) = ^ m + a)+.(n Cm -a) ^ m=12> (34) 

This condition can be intuitively understood by the following picture: The linear chain 
is driven by the washboard wave which has "momentum" (i.e., wave number) 27rc and 
"energy" (i.e., frequency) 2nv. If the driving is sufficiently strong, phonons with wave 



number k± = ircm±q are excited parametrically by m washboard waves. The condition (pl|) 
follows from the conservation laws of momentum and energy. It depends on a continuous 
parameter q contrary to the resonance condition ([3(]). Of course the periodic boundary 
condition (0) restrict the value of q to N different values corresponding to N pairs of phonons. 
For c = 1 the driving is uniform and ( |33] ) becomes identical to ( |33"D where c^o is the frequency 
of the pair of phonons with the same wave number but propagating in opposite directions 
and forming a standing wave. For weak damping the threshold for the m-th order parametric 
resonance is proportional to j 1 /" 1 (see in appendix ||) . 

The maxima and minima of v^(q) define resonance intervals which successively appear for 
increasing b starting with the first-order parametric resonance. There is a relation between 
these instability intervals and the superharmonic resonances: 

< > «£• (35) 

Thus the first-order parametric resonance occurs for sliding velocities larger than the velocity 
vf of the main resonance. For k — > the minimum of is given by v^. It is possible that 



a superharmonic resonance peak sits in a parametric resonance interval, i.e., mmv^ < v f < 
maxf^ for n ^ m > 1. Also overlaps of instability intervals are possible. 

For sliding velocities in an instability interval, the amplitudes of two phonons which are 
selected by v = v ^ increases exponentially. This exponential increase is saturated by the 
nonlinear terms omitted in flTT|). Their effects will be discussed in the following subsections. 
The increase and saturation of the amplitudes also occur in Fk- We expect that parametric 
resonance does not lead to a peak in Fk{v) but to something which looks like a more or 
less flat table mountain. Furthermore, we expect quasiperiodic motion because the phonon 
frequencies uj(tc cm + q) and u(7rcm — q) are in general incommensurate except for c = 1. 



B. Periodic and quasiperiodic solutions 



From section |11 A| we know that the nonlinear equation of motion (^) is invariant under 



discrete translation in time. Solutions with the same symmetry should therefore exist. They 
can be described by a one-parameter function called dynamical hull function g£>\ 

= 9D(vt + cj), with g D (x + 1) = g D (x). (36) 

The dynamical hull function is very similar to the static hull function g$ introduced in 
Ref. 0. In fact for v — > the dynamical hull function becomes identical to the static one. 
It is a solution of 

v 2 9d(x) + lvg' D (x) = g D (x - c) + g D (x + c) - 

-{2 + K)g D {x) + 

+6sin27r[x + g D (x)]. (37) 

With the Fourier ansatz 

oo 

g D [x)= £ A m e 2mmx , A_ m = A* m . (38) 

m=— oo 

this differential-delay equation turns into a set of nonlinear algebraic equations: 

L m A m = b J 1 sin j 2vrx + 2tt £ Ae 2 ™ 1 * J e~ 27Timx dx. (39) 



(=— oo 



We solve this equation numerically by truncating the number of Fourier coefficients. Fig- 
ure |](b) shows the result for an increasing sequence of values of b. 

We see that the dynamical hull function is not uniquely defined if b exceeds some thresh- 
old. The resonance peaks starting with the main resonance develop loops leading to bista- 
bility and hystereses. This is caused by the fact that the oscillation amplitude of the particle 
becomes of the same order as the periodicity of the external potential. These loops are also 
found in a single Tomlinson oscillator |10| , |TT| . 

The nonlinear terms lead also to a reduction of the height of the resonance peaks com- 
pared to what is expected from the linear theory [see e.g. Fig. [|(a)]. Because the peak 
height in the linear theory is proportional to the inverse damping constant [see fl23|)1, we 
interpret this result as an increase of the effective damping constant. The reason for this 
increase is the opening of additional channels of energy dissipation due to nonlinear phonon 
coupling. 

In order to investigate the stability of the solution described by the dynamical hull 
function, one has to study the equation of motion linearized around (PE|). That is, we 
consider only those terms which are linear in 8£j = £j — go{vt + cj). We get 

<k, + = K, ■ : + <K, ■ : - (2 + /.'Ms, • 

+2tt6 cos(2tt[^ + cj + g D (vt + cj)]) 5£j. (40) 

With the help of the well-known Floquet-Bloch ansatz for 8£j (see appendix [B]) we solve this 
equation numerically with truncated Fourier expansions of go and 8£j. 

There are two types of instabilities of the hull function. The first type is caused by 
the bistability of the resonance loop because there is always a separating unstable solution 
between the stable ones (called a saddle in the terminology of bifurcation theory). The 
other instability type is caused by parametric resonance already discussed for the linearized 
equation of motion ( pf ) in Sec. [IV A3| . From the stability analysis in appendix |B| we know 



that the most destabilizing disturbance contains two generically incommensurable frequen- 
cies (2ttv and Q). This is confirmed numerically for flip]) . At these instabilities the system 
selects a state which leads most often to an enhanced kinetic friction. An example for this 
enhancement is shown in Fig. §(c). The most prominent enhancement occurs for first-order 



parametric resonance. 

Figure ||(a) shows the quasiperiodic motion expected from Sec. [IV A 3| . It shows a stro- 



boscope map of the lateral force in time steps given the ratio of the lattice constant of the 
external potential divided by the velocity. The lateral force is given by 

K N 
3=1 

The dynamical hull function corresponds to a fixed point of the map 

fn+l = T(f n ). 

Quasiperiodic motion leads to a closed one- dimensional attractor of the map. For c = 1 the 
motion is periodic leading to a limit cycle of period two. 

C. Irregular behavior 

For increasing driving strength b the quasiperiodic solutions caused by parametric reso- 
nance becomes unstable and the motion become more and more irregular and chaotic [see 
Fig. |5](b)]. It is presumably spatio-temporally chaotic as in the Frenkel-Kontorova model 
| 36fl . That is, the number of positive Lyapunov exponents is an extensive quantity |37 



i.e., 



it is proportional to N. We call this strongly mixing state a fluid- sliding state even though 
the harmonic interaction in the model excludes on average changes in the particle order. 
Nevertheless the state is fluid-like because the amplitudes of the particle motion are of the 
same order as the lattice constant c. 

The kinetic friction F% from simulations are shown in Figure ^. In each simulation the 
sliding velocity was increased in tiny steps from zero up to a certain value and then decreased 
back to zero. Beyond first-order parametric resonance (i.e., v > maxff) large hysteresis 
loops occur leading to bistability between fluid-sliding states and the solid-sliding state. 

Characteristic peaks appear in Fk(v) in the regime of fluid-sliding states beyond the 
main resonance. What is the nature of these peaks? Looking at the spatially and temporally 
resolved dynamics we were not able to see any difference between the dynamics inside and 
outside the peaks. The motion is strongly irregular. The only exceptions are the large peaks 



at v = 1.06 and v = 0.95 for c = 144/233 and c = 1, respectively. Here coherent motion 
appears which is investigated in the following subsection. 

The spatio-temporal Fourier transformation of £j(t) reveals a structured distribution in 
the (wave-number, frequency) space which is mostly concentrated on the curve defined by 
the dispersion relation flTSp. Therefore the strongly chaotic motion in the fluid-sliding state 
is mainly caused by weak nonlinear interactions between phonons. 

Here we will give only some qualitative arguments for the occurrence of the peaks beyond 
the main resonance. The argument is based on the picture of multi-phonon scattering 



processes already used in subsection |IV A| to explain parametric resonance. The most general 
process which drives the phonon system is the decay of m washboard waves into n different 
phonons. From energy and momentum conservation we get 



and 



(41a) 



27r(cm + u) = kj , with u integer. (41b) 
j'=i 

The integer u denotes possible Umklapp processes. The strength of the decay processes 
scales like b m . The most important processes are those with small n and m. Superharmonic 
resonance (i.e., n — 1) and parametric resonance (i.e., n — 2) are the linear cases. For n > 2 
the decay processes are nonlinear. Beyond first-order parametric resonance (i.e., v > maxvf) 
the solid-sliding state is therefore insensitive to these processes. They become active only 
if at least some phonons are excited. Then they excite additional phonons leading, to a 
self-sustained process that converts macroscopic, uniform sliding into microscopic, irregular 
motion. 

An important quantity characterizing the efficiency of phonon driving due to washboard 
waves is the density of combinations of wave vectors fulfilling fl41| ) for sliding velocities from 
the interval of [v,v + dv\. We call this density the density of decay channels. We expect 
that characteristic structures of this density, such as peaks, also appear in Fk(v). 

Because of the constant wave-number density due to periodic boundary conditions, the 
density of decay channels of an n-wave decay process of order m reads 
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^ ' 

where 5 is Dirac's delta function and k n in is substituted due to ( |41b|) . 
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Extrema of lead to singularities in p 7 ^ which are generalizations of van-Hove singular- 
ities well-known from phonon spectra ||38|| . In appendix |Cj we have calculated the extrema. 
For the parameters of Fig. ^| the most important extrema are tabulated in table |. The 
scaling behavior of pj^ near an extremum depends on its type. The most prominent singu- 
larities are the diverging ones. There are two cases of diverging singularities. The first one 
occurs only for n = 2 where the density scales like \v — f c | -1 ^ 2 near maxima and minima of 
v^. For n > 2 logarithmic singularities occur near certain types of saddles. A saddle can be 
characterized by the numbers n+ and n_ which count the linearly independent directions in 
which increases or decreases, respectively. Logarithmic singularities occur only if both 
numbers are odd (see appendix [D]). 

For the parameters of Fig. ^| we have calculated numerically the density of the most 
important decay processes. The results are shown in figure [7[ Together with table | we are 
able to identify all peaks in Fk{v) that occur beyond the main resonance. 

Concerning the strength of friction this is only a qualitative analysis. The reason for that 
is twofold. First, the decay channels are not equally fed and, second, we do not know the 
effective damping constant of the phonon modes due to nonlinear phonon scattering. This 
information would be needed in order to obtain a quantitative theory of F K {v). 



D. Coherent motions and solitons 

As already mentioned in the previous subsection, the motion is not entirely chaotic near 
the large peaks at v = 1.06 and v = 0.95 for c = 144/233 and c = 1, respectively. Coherent 
waves occur, as can be seen in Figure |8|. A spatio-temporal Fourier analysis shows a peak 
which does not sit on the dispersion relation but near the point which is determined by the 
wave number and frequency of the phonon of the decay process that belongs to the absolute 
maximum of vf. 

In the example shown in Fig. || we observe a modulation of the amplitude of the fast 
oscillation. We have extracted the amplitude of this envelope from the data of the simula- 
tion. Figure ^| shows the spatio-temporal evolution of this amplitude. The regular motion of 
the equidistant maxima of the envelope in Fig. || appears as parallel white stripes in Fig. |9|. 



On this background pattern three dark stripes are visible. They are a kind of dark solitons 
well-known from the nonlinear Schrodinger equation. Since they propagate with different 
velocities, collisions occur which are nondestructive like in the nonlinear Schrodinger equa- 
tion. Thus the FKT model seems to behave like a conservative, undamped and undriven 
system. But the waves emitted by the solitons are strongly damped. 



V. CONCLUSION 

We have studied the dynamical properties of the FKT model driven at a constant sliding 
velocity v. The model describes in a rudimentary way the interface between two atomically 
flat sliding surfaces. An advantage of this model is that some basic questions concerning 
dry friction can be studied analytically. 

Since dry friction operates mostly far from equilibrium, we are in principle not able to 
calculate macroscopic properties by using straightforward recipes. In the sliding process of 
two solids, the uniform motion is turned into microscopic, irregular motion called heat. The 
rate of dissipation determines the kinetic friction Fx- What are the mechanisms of this 
transfer from macroscopic to microscopic motion? How effective are these mechanisms? 

The two basic dissipation channels in the FKT model are the damping terms which are 
proportional to the velocities of the interface particles relative to the sliding bodies. The 
basic mechanisms of energy transfer are multi-phonon decay processes where m so-called 
washboard waves excite n phonons of the chain of interface particles. The washboard is the 
spatially periodic potential describing the atomically flat surface of the lower body. It acts 
as a wave (wave length = surface lattice constant, frequency = sliding velocity divided by 
wave length) in the frame of the upper body to which the surface particles are connected 
by springs. The washboard is driving the surface particles. For small driving strength b and 
large damping the process with m = 1 and n = 1 is the most dominant one. It simply leads 
to a resonance peak in the kinetic friction Fk(v) given by A broad peak is exactly 



what is found in boundary lubrication [33]. For large velocities the washboard shakes the 
particles so rapidly that its influence is averaged out. Therefore the interaction becomes 
effectively zero. The interface motion freezes out forming the solid-sliding state. In this 



state the kinetic friction is proportional to the sliding velocity. Note that in Figs. ^ and |6| 
this part of the kinetic friction is always omitted. 

For increasing b or decreasing damping superharmonic resonance peaks occur. They 
are caused by processes with m > 1 and n — 1. The superharmonic resonance peaks are 
smaller than the main resonance peak and occur for velocities below the main resonance. 
As long as the oscillation amplitudes of the interface particles are much less than the wash- 
board wave-length, the main resonance and superharmonic resonances can by approximated 
analytically very well from the linearized equation of motion. Otherwise processes with 
n > 1 become important leading to nonlinear phonon damping. This increases the effective 
damping constant resulting in a decrease of the resonance peaks in accordance with (|23|). 

The decay process with m = 1 and n = 2 is responsible for an instability (first-order 
parametric resonance) of the solid-sliding state if b exceeds some threshold which is pro- 
portional to the damping constant. This instability occurs for a whole interval of v beyond 
main resonance. The amplitudes of both phonons increase. Nonlinearities in the equation 
of motion eventually saturate the growth due to nonlinear phonon damping. This leads to 
quasiperiodic motion because the phonon frequencies are in general incommensurate. The 
kinetic friction in the parametric resonance window can be larger than at the main resonance. 

A further increase of b leads to bistability between the solid-sliding state and the so-called 
fluid-sliding state which is characterized by spatio-temporal chaos. The kinetic friction of 
this state exhibits some peaks which can be qualitatively understood by the density of 
decay channels for processes with m = 1 and n = 2, 3, and 4. In the regime of fluid-sliding 
states a small velocity interval is embedded where coherent motion including nondestructive 
collisions of dark envelope solitons occurs. 

It should be noted that the behavior beyond main resonance does not depend on the 
commensurability of the ratio c of the surface lattice constants of the sliding bodies, contrary 
to what is found for small sliding velocities. 

The behavior for small sliding velocities is dominated by stick-slip motion of the interface 
particles. In the overdamped limit every Q-th particle (Q is the denominator of c) jumps 
into the neighboring potential well leading to F$ — Fk{v — > 0) ~ 1/Q- Thus the kinetic 
friction Fk(v — > 0) is identical with the static friction Fg in the incommensurate case. The 



kinetic friction is proportional to the sliding velocity if b < bf. In any case Fk{v) is an 
increasing function at v — 0. For finite damping constants particles may jump to more 
distant potential wells leading to kinetic friction which is less than the kinetic friction in the 
overdamped limit. Furthermore, this stick-slip motion becomes more irregular. 

In order to compare the dynamical behavior of the FKT model with the BK model we 
have to note that most studies at the BK model are done for k < 1 whereas in our study we 
have always chosen k = 1. The behavior of the BK model is different for k = 1 compared to 
k C 1. For example, the solitonic motion of a group of a few blocks found by Schmittbuhl 



et al. |20| disappears if k is too large WA - F° r K — 1-2 Sarkadei and Jacobs p3| have 
found regular periodic waves. They travel from the open end into the chain and annihilate 
somewhere in the middle. But contrary to the FKT model, this behavior does not lead to 
resonance peaks in the kinetic friction ||22,p3|. The reason for that may be the fact that in 
the BK model energy is dissipated during the sticking phases of the blocks and not during 
sliding. 

There are several theoretical studies where the kinetic friction is investi- 

gated in a system which is driven by a constant force and not by a constant velocity as in 
our study of the FKT model. The averaged velocity of the center of mass is calculated. All 
these studies show fluid-sliding states with enhanced friction and solid-sliding states where 
the kinetic friction is given by the sliding velocity times the damping constant. Hysteresis 
between the solid-sliding state and fluid-sliding state occurs only if the temperature is small 



enough |[42|| . Therefore we expect that the strong hysteresis in the FKT model disappears 
for finite and large enough temperatures. 

A model which goes beyond the FKT model has been studied by Sokoloff P, |29| , f43 



Instead of one layer it has several layers coupled harmonically. Thus the driving due to the 
washboard wave is no longer uniform. Sokoloff has calculated the kinetic friction analytically 
for a prescribed periodic motion of the layer which feels the periodic potential. The result 
is therefore similar our results for the main resonance and superharmonic resonance (if 
the prescribed motion has higher harmonics). This kind of approach is not able to deal 
with instabilities like parametric resonances or other multi-phonon processes. It would be 
worthwhile to study such an extension of the FKT model from our point of view. 



In a forthcoming paper we will study macroscopic stick-slip motion of the FKT model 
by coupling the upper sliding body via a spring at a support which moves with constant 
velocity. We are interested in the following questions. What is the interplay between static 
friction and kinetic friction? What happens during the transitions from stick to slip and vice 
versa? Which (meta)stable states does the system select after the slip-to-stick transition? 
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APPENDIX A: QUASISTATIC SLIDING IN THE TOMLINSON LIMIT AND 

FOR c = 1 

In the Tomlinson limit (i.e., n,b — > oo) and for c = 1 the static hull function g$ can be 
given in parametric form because the next-neighbor coupling in @ can be dropped: || 



In order to calculate Fjc in the quasistatic limit v — ► we have to distinguish between 
b<b* = k/(2tt) and b > bff . 

In the first case we can apply ( |T0D with 7 c g defined by (|TT|). To do this the integral 
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9s(r) = -sin27rr, x(r) 



r-gs{r). 



(Al) 



Jo (dgs/dx) 



2 dx has to be calculated. Using ( |A~1| ) we get 




(A2) 



Thus 



/ 



1 



(A3) 



7eff = 1L + 7 




1 



For b > b£ the hull function g$ is no longer a unique function of x. There is a saddle-node 
bifurcation at the relative maximum of x(t) at 

he 

2,tyt§]\[ = arccos — — — . (A4) 
2i\b 

For the same value of x there is at least a second solution defining a r 2 implicitly by 

t% - #s(t 2 ) = t sn - g s (T SN ). (A5) 
This equation can be solved perturbatively near b = bf and n/b — > leading to 



+^f-(^-^)V /2 + W-^) 7/2 ) (A6) 



350 \k 



and 



7T 



/2nK k 1 /n/t\ 3 / 2 
2. T2 = 2m ---^— + — --(-) + 

+ ^(¥) 2 + ° ((K/,,)5/2) - (A7) 
with n = 1,2, . . ., respectively. 

Because x corresponds to the relative position of the sliding bodies, an increase of x 
beyond the relative maximum yields a jump of one particle into a state corresponding to the 
solution of (|A5|) . In order to calculate the kinetic friction it is most convenient to investigate 
the difference of the potential energy before and after the jump: 

b 

AV = — [cos27rTsjv — cos27tt 2 ] + 
2% 

+ ^[9s(tsn) 2 - 9s(r 2 )}. (A8) 



This expression can be simplified by using ( |Alf ) and (TA4 ): 



k / 2nb 



AV = —r 1 COS 27TT 2 . (A9) 



6ir 2 V K 



In accordance with (^) the kinetic friction reads 

AV 

F K (v^0) = —N, (A10) 

where (n) is the number of potential wells a particle jumps on average. For b —>■ bf a particle 
can jump only into the next potential well. Thus (n) = 1 and we get with (|A6[) 



JM 2k dk z 

2UllT \b-b?r + o((b-b«n (ah) 



175k 3 

In the limit n/b — > we get with (|A7|) 



F K (v - 0) ( (n)K 



N \ 2b 

+ ^r{b) + °^i h ))\- ( A12 ) 

These equations generalize the results of Ref. [ITOf. 



APPENDIX B: PERTURBATION THEORY FOR THE LINEARIZED 

EQUATION OF MOTION 



The linearized equation of motion fli~7|) can be treated analytically in the framework of 



a perturbation theory with b as the smallness parameter. 

First the solution is calculated in the following way. We start with (^7]) and expand the 
solution: 

oo 

An = EAn/- (Bl) 
1=1 

Inserting this ansatz into ( p7[ ) leads to a recursion equation for the coefficients A m /. 

i 

A m ,l = 7TT (<5m,-l — 5m,l) (B2) 

2Llr, 



and 



7T 

A m ,i+i = t — (Am+ij + A m -i } i). (B3) 

J-lm. 



The solution of this recursion scheme is 



ITT 1 ' 1 1 1 1 



Thus the only nonzero coefficients are 

A m , m +2n, with m + n / 0, m,n = 0,1, ... (B5) 



Therefore the leading term of the Fourier component A m is of order b n 



- x 

i=l -^v 

2 m f_Z,\2 



Using (|29| ) we get for the kinetic friction 

oo I l—m 

-S = ^ 2 vj £ b21 E m 'E A m , m+2n ^ i2i _ m _ 2n , (B7) 

iV 1=1 m=l n=0 

which is an even function in b as expected from symmetry. 

In order to study the stability of the solutions of we make the following ansatz in 
accordance with the well-known Floquet-Bloch theorem : 

oo 

= E a m e 2nim{vt+cj)+iqj+xt , (B8) 

m=— oo 

where g is an arbitrary given wave number (with qN an integer multiple of 2ir) and A a 
complex number which is an eigenvalue of: 

D m a m - 7r6(a m „i + a m+1 ) = 0, (B9) 

with 

D m = (2-Kivrn + A) 2 + 7(27m;m + A) + u 2 (2ircm + q). (BIO) 



Note that ( |B9|) is invariant under the transformation A — > \ + 2niv, q — > q + 2nc, a m — > a m+1 . 
A solution of (|T7|) is stable as long as the real parts of all eigenvalues are negative. In order 
to calculate the instability threshold b c it is more convenient to solve ( p9|) for A = ifl where 
b and Q are the unknown variables. The diagonal elements turn into 

D rn = R m + i 7 I m , (Bll) 

with 

R m = uj 2 (2ircm + q) - I 2 m , I m = 2-Kvm + H (B12) 

The real solution b of ( |B"9| ) with the smallest absolute value defines a function of q. Its 
absolute minimum gives the instability threshold b c . It is always zero in the undamped case 
(i.e., 7 = 0) because at least one D m is zero. In the case of damping all -D m 's are nonzero. 



Parametric resonance of order n corresponds to the fact that two diagonal elements, say 
D and D n , go to zero for b — > and 7^0. In order to calculate the instability threshold 
we assume that a$ and a n are given. Next we express a±\ and a n ±\ in terms of ao and a n : 
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(B13) 
(B14) 
(B15) 
(B16) 



where 



a"; 



7r6 D n+ i 



tt6 D m 

Eq. ( P39|) for m = and m = n turns into a set of linear equations for ao and a n . It has 



(B17) 



nontrivial solutions if the determinant 

n d~ 2 d n ~ l 
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is zero. The assumption 



R ,R n = O(b 2 ), R m = O(b ), m^0,n, 



yields 



'n-l 



HD j + 0(b 2 )) -^I I n + O(b 2 ) =(7rb) 2n . 



\m=l 



Thus in leading order of 7 we obtain 



/ n-l \ V" 

= 7\/-Wn II l-^il 
V m=l / 



(B19) 



(B20) 



(B21) 



The product term does not appear for n = 1. Bearing in mind that IqI u has to be negative 
and u(k) is an even function we immediately get the parametric resonance condition (31) 
by solving R = R n = 0. 



APPENDIX C: EXTREMA OF i£ 



Before calculating the extrema of t>™ defined by PT|) we eliminate k n by using ( 41b|) . 
Thus is a function of k±, . . . , k n -\. An extremum is given by 

27rm9 fc .< = u/(%) - u/(A:„) = 0, j = 1, . . . , n - 1, (CI) 

where 

eL; sin A; 

" (fc) = ^ = ^(fc) (C2) 
is the derivative of the dispersion relation (|T8|). For an extremum all u'(kj) have to be 
the same. Because u'(k) is a periodic function, there exist two solutions and fc+ in the 
interval [— it, 7r). For a given value of we can calculate the corresponding by solving 
uj'{kJ) = u/(k+): 

t • ^ / 2-(/c + 2) cosfc + ^ 

= s lg n (it,) arccos ^ +2 _ 2 . , (G3) 



where sign is the signum function. Therefore an extremum is given by 

puj(k_) + (n -p)w(fc+) 



(C4) 



2vrm 

where p and n — p are the numbers of wave vectors which have the wave numbers and 
respectively. 

The wave number fc_ is not arbitrary because momentum conservation ( [41b|) has to be 
satisfied: 

pk^ + (n — p)k + = + u), with it integer. (C5) 



Eliminating in (|U5|) by using ( |C3|) leads to a nonlinear algebraic equation for k + which 
can be solved only numerically except for p = 0. In this case all wave vectors are identical. 
That is, 

kj — k + — 2n — — , j — 1, . . . , n. (C6) 

The type of the extremum is determined by the eigenvalues of 

2ixmd k] d h v n m = u"(kj)8u + u"(kn), (C7) 



where u"(k) is the second derivative of the dispersion relation (0). For p = the eigenvalues 
are 

Ai = nuj"(k + ), Xj = u>"(k+), j = 2, . . . , n — 1. (C8) 

Thus the extremum is either a maximum or minimum depending on the sign of uj" . In the 
general case we have two eigenvalues given by 



^ = 2 ' ( C9 ) 

where 

u'l=u}"(k±), (CIO) 
p — 1 times the eigenvalue A = lo"__ , and n — p — 2 times the eigenvalue A = ui'L . 



APPENDIX D: SCALING OF p™ NEAR EXTREMA OF v. 



The scaling behavior of p 7 ^ near an extremum of v^{kx, . . . , k n -x) can be easily ob- 
tained after a linear transformation of the coordinate system k\, . . . , k n _\ into a new system 
xx, ■ ■ ■ i sc n --i which is given by the eigendirections of the matrix defined by the second deriva- 
tives of calculated at the extremum. Up to second order we get 

*C = «c + ( D1 ) 

3=1 3=1 

where n + (n_) is the number of positive (negative) eigenvalues of that matrix. A non-singular 
matrix leads to n + + n_ = n — 1, otherwise n + + n_ < n — 1. 

The scaling behavior near an extremum is determined by the states nearby. Thus we 
perform the integration (|42| ) only in a spherical domain of radius R: 

Pm ~ /* - ^ - r + + rlyl^rl^dr+dr., (D2) 
where r± = Ej=i ^f- O ne integration yields 

ft; - u c - r 2 ) n+/2 - 1 r"-~ 1 rfr_. (D3) 



m 







In the case of either a relative minimum or maximum (i.e., n_ = or n + = 0) we get 
therefore 



pi ~ (v - v c ) n + Z 2 " 1 or ^( Vc _ v )"-/2-i. (D4) 

In the case of a saddle (i.e., n± 7^ 0) we perform the second integration ( P3|) by using the 
Euler substitution 



t = r_ + y z + ri, with z = v — v c , 

which yields 

p n m ~ J {t 2 +z) n+ ~ 1 {t 2 -z) n — H l - n +- n ~dt. (D5) 
Neglecting the contribution from the upper limit we get 



j+=o j_=o V J+ / \ J- 



x 



X2— l I(n+-7i_ -2-2j+-2j_), (D6) 



where 



. lnz, s = 

7(5) = <| . (D7) 

1, s^O 



The logarithmic term appears only if 



n± — 1 



2(j + +j_+l) = n + +n_ 

is not zero. It is easy to show that C = — (— l) n± C holds. Thus 



11 _]_ + n _ 



, In If — v c \, n± are odd 
P n m ~\v-v c \-^-H . (D9) 



1 , otherwise 
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FIGURES 



FIG. 1. A schematical sketch of the Frenkel-Kontorova-Tomlinson (FKT) model. 

FIG. 2. The friction force in the limit of quasistatic sliding (i.e., v — ► 0) for k = 1. The solid 
lines show the kinetic friction Fk in the overdamped limit (i.e., 7 — > 00). The dashed lines denote 
the static friction i^. F™ ax = frjV is the overall upper bound of static friction. For c = 1 the 
analytically obtained approximations of Fk, ( All| ) and (|l(]), are shown by dotted lines. 



FIG. 3. The kinetic friction Fk for small sliding velocities. Parameters are k = 1, N = 233, 
M = 144, 6 = 2/vr w 0.637. Using definition (j^) the data has been obtained by numerical 
integration of (Hj). Each data point is the average over 500/u time units. In order to avoid averaging 
during transients the system has been evolved over 10 3 time units before starting averaging. The 
last state has been chosen as the initial value for the next data point with increased v. In order to 
get some information about the low-frequency fluctuations the averaging interval is chunked into 
ten pieces, each time units long, yielding ten values of Fk- The standard deviation of them 
are denoted by error bars. 

FIG. 4. Main resonance, superharmonic and parametric resonance for k = 1, N = 233, and 
M = 144. Resonance velocities of main resonance and superharmonic resonances in accordance 
with (|30|): vf ~ 0.337, ~ 0.134, v% ~ 0.071, vf ~ 0.089. Parametric resonance intervals in 
accordance with (|34|): First order: [0.480,0.614]. Second order: [0.197,0.337]. Part (a) shows a 
comparison between the approximation (^) (dotted line) and the numerically obtained solution of 
equation of the dynamical hull function (|37| ) (solid line) for 7 = 0.5 and b = 1/tt ~ 0.318. Part (b) 
shows different solutions of ([}?]) for 7 = 0.1 and 2irb = 0.1, 0.2, 0.5, 1 and 2. Solid (dotted) lines 
indicate stable (unstable) solutions. Part (c) compares the hull function solution (dashed line) 
with the full simulation of (|5|) (solid line and squares for increasing v, dotted line and triangles for 
decreasing v) for 7 = 0.1 and b = l/(2-/r) ~ 0.159. Details of the simulation are the same as in 
Fig.§ 



FIG. 5. Stroboscopic map of quasiper iodic and chaotic motion caused by first-order paramet- 
ric resonance. The parameters are k = 1, N = 233, M = 144, j = 0.1, v = 0.5, and (a) 
b = l/(2vr) 0.159, (b) b = 3/(2vr) 0.478. /„ is the lateral force -k/NJ2j £j at time t = n/v. 

FIG. 6. The kinetic friction F K (v). Parameters are k = 1, N = 233, b = 3/(2vr) 0.478, 
7 = 0.1, and (a) c = 144/233, (b) c = 1. Squares (triangles) connected by solid (dashed) lines 
denote the numerically obtained values for increasing (decreasing) sliding velocity v. The averaging 
procedure is the same as described in the caption of Fig. |3|. 

FIG. 7. Densities of the most important decay channels beyond the main resonance. The 
parameters are k = 1 and (a) c = (\/5 -l)/2w 0.618, (b) c = 1. 

FIG. 8. Coherent motion. The first five particles are shown. The parameters are k = 1, 
N = 233, M = 144, 7 = 0.1, b = 3/(2vr) » 0.478, and u = 1.06. 

FIG. 9. Dark envelope solitons embedded into the coherent motion seen in Fig. ||. The param- 
eters are the same as in Fig. |8| The grey scale from dark to bright denotes the strength of the 
amplitude of the fast oscillations of the particles. The coherent motion is indicated by the regular 
stripes. There are three dark solitons traveling with three different velocities. At j ~ 100 and 
t « 190 two of them interact nondestructively. 
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Fig. 7(a) — M. Weiss and F.J. Elmer 




Fig. 7(b) — M. Weiss and F.J. Elmer 
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Fig. 4(c) — M. Weiss and F.J. Elmer 
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